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Abstract 

We present an exact solution for the time-dependent Stokes problem of 
an infinite cylinder of radius r = a in a fluid with harmonic boundary 
conditions at infinity. This is a 3-dimensionaI problem but, because of 
translational invariance along the axis of the cylinder it effectively reduces 
to a 2-dimensional one. The Stokes problem being a linear reduction of the 
full Navier-Stokes equations, we show how to satisfy the no-slip boundary 
condition at the cylinder surface and the harmonic boundary condition 
at infinity, exhibit the full velocity field for radius r > a, and discuss the 
nature of the solutions for the specific case of air at sea level. 



1 Introduction: Navier-Stokes Equations 



We denote a position in space by x and time by t. We describe an arbitrary fluid 
flow by its local velocity v{x, t), density g{x, t), and pressure p{x, t) and denote the 
total derivative by 

Dg{x,t) dg{x,t) A dg{x,t) 



dt 



i=l 



dxi 



For the sake of clarity, we will omit the position and time dependency and simply 
write Q = g{x, t). In an incompressible fluid of dynamic vis cosity 1^ = (j,/ g and with 
volume force F the Navier-Stokes equations are given by (jPantonl . l2005l ) 



Dv 



dv 



dt go 
while incompressibility is explicitly taken care of by 

div V = . 



gradp + Av 



Qo 



(2) 



(3) 



Finally, we will use the no-slip boundary condition, meaning that the velocity van- 
ishes, i.e., V = 0, at the surface of an inflnite cylinder whose axis we take parallel 
to the 2;-axis. The cylinder radius equals r = a. 



constant 


value 


unit 




15.11 ■ 10-^ 


m^/s 


Qo 


1.204 


kg/m^ 



Table 1: Measured values for the relevant material dependent constants for dry air 
at temperature 20°C and pressure po = 1013 hPa. 



The problem we are going to solve exactly is that of inflnite cylinder in the above 
described fluid with harmonic boundary conditions at inflnity. Since the cylinder 
itself is inflnitely long we can simplify the problem to a 2-dimensional one. This we 
assume throughout what follows. 



Solution of the Stokes Equations for an Infinite 
Cylinder in a Viscous Harmonic Flow 



Stoked (Il85ll ) has already calculated the effect of internal friction of a harmonically 



moving fluid with angular frequency u on the motion of pendulums with radius a 
in an approximation that is valid for a^fujJvQ ^ 1. In the following we will give an 
analytical solution for the velocity fleld that is valid for all r > a and a range of 
frequencies uj with Reynolds number Re -C 1. 
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In practical work the stream velocity may well vary arbitrarily in time but, to 
simplify the problem, we solve it for a harmonic flow field. Hence the stream velocity 
at infinite distance from the cylinder is 




X — —CO X = oo 



Figure 1: Cross-section through a cylinder with radius a in a flow field that is 
harmonic at infinite distance from the cylinder. Angle 6 and distance r indicate the 
polar coordinates used here. 



For Re <^ 1 the non-linear term {v ■ 'V)v is negligible and we can use the time- 
dependent Stokes equation 

dv gradp 

— = + uoAv. (5) 

at go 

It is linear in v. Constant volume forces such as gravity do not play any role here 
since they only cause additional gravitational pressure and do not change the form 
of the equations as they can be incorporated into gradp. 



2.1 Boundary Conditions 

The boundary condition (jl]) at infinite distance from the cylinder is consistent with 
to the time-dependent Stokes equation ([5]) as well as with the Navier-Stokes equa- 
tions ([2]). This is important because at some distance from the cylinder the time- 
dependent Stokes equation is not valid any more. 
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On the cylinder surface we have no-shp boundary conditions, 

v{a,t) = 0. (6) 

Using the boundary condition (jl]) at infinity we can calculate the pressure at infinite 
distance from the origin. As the velocity at infinity is homogeneous 

Av^ = (7) 

and {v ■ 'V)v has been dropped we get, using ([2]), 

dvr 



gradpoo = -^'0- 



dt 

In view of Fig. [1] the boundary condition for the pressure at infinity is therefore 

Poo = ixQouvoe'''^* . (9) 

2.2 Solution of the Stokes Equation 

Applying the divergence to both sides of the time-dependent Stokes equation 
and taking advantage of we get 

ddivv div gradp 

— - — = \-uAdivv (10) 

dt go 

and as div v = we therefore find for the pressure 

Ap = 0. (11) 

We first give a general solution to ( fTTj) that agrees with the boundary condition ([9]) 
and the symmetry of the problem. 

In cylindrical coordinates, the boundary condition iQ for the pressure at infinity is 

p(oo, 9) = ir cos OgoUVoe''"'^^ (12) 

and ffTTl) becomes 

Ap = ^ + + — ^ + ^ = (13) 

We give a general solution using the ansatz 

p{r,e,t) = R{r)P{e)T{t) . (14) 

We get from f|T3l) . as R{r) and P{0) are functions depending on different variables, 
with a constant k independent of r and 9, 

^ + = (15) 
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and 



d'^R r dR , , 
R dr^ R dr 



0. 



The general solution for R{r) of (fT6l) is 



R{r) = ^-A,,r^ 

while for P{d) 

P(e) = Bi^k cos{ke) + B2,k sm{ke) . 
The pressure must be continuous and hence 

P{e) = P{e + 27r) 
so that /c = 0,l,2,....Asa general solution to ( fT3l) we get 



(16) 

(17) 

(18) 
(19) 



k=0 



[Bi,k cos (ke) + 52,fc sin (A;^)] T{t) . (20) 



To comply with f[T^ the solution cannot contain any power of r'^ with k > 1 since 
the functions cos{k6) and sin(fc^) are linearly independent for all /c = 0, 1, 2, ... . As 
the Stokes equation is linear we set T{t) = exp{—iujt). Expressions with k = are 
irrelevant since they only result in an additive constant to the pressure. Possible 
solutions are therefore 



p(r, = I'^-j [Ci,k cos (kO) + C2,k sin {kO)] + irQ^uvQ cos e 



. k=l 



The problem is symmetric with respect to the x-axis and therefore 



(21) 



(22) 



so that solutions depending on sin(A;^) must be zero. As we get another solution of 
the time-dependent Stokes equations by applying the transformation v —v and 
p — >■ —p, which in our case is equivalent to mirroring with respect to the y-axis, it 
holds that p{—x,y) = —p{x,y) or, in polar coordinates, 



p{r,e,t) = -p{r,n - 9,t) . 



Since 



cos[A; (tt - 9)] 



cos{k6) for k even 
- cos(A;^^) for k odd. 



(23) 



(24) 



the solutions depending on cos{k6) with even k must vanish. The remaining most 
general solution that complies with the boundary conditions is 



p{r,9,t) 



oo ^ 

Y,—CkCos{ke) + 



irQQOJVQ cos 9 



k=l 



-iujt 



(25) 
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We do not know yet the boundary condition for the pressure at the surface of the 
cyhnder as it results from the flow field. To simplify further calculations we use the 
simple ansatz 



p(r,9,t) 



C 



cos ( 



-iujt 



(26) 



and verify that it is consistent with the Stokes equations 



Since we are in two dimensions, the incompressibility condition reads in polar coor- 
dinates 



1 d irvr 

awv = — 

r or 



r do 







and the time-dependent Stokes equations are given by 



dvr 
'dt 

dve 
dt 



+ V 



1 9p 
Qq dr 
1 dp 
tqq do 



^ d'^Vr 1 d'^Vr 1 dVr 2 dVf) 



+ 



+ 



d'^ve I d'^ve 1 dve 2 dvr 

Qj,2 y,2 QQ2 Qy, y,2 QQ 



(27) 

(28) 
(29) 



We now substitute the pressure ansatz (!26|) . With the help of (!27|) the above equa- 
tions decouple so as to give an equation for Vr alone, 



dvr 
'dt 



C 
gor^ 



tUJVQ 



cost 



—iuit 



+ V 



d'^Vr 
Qy,2 



1 d'^Vr 3 dVr 
+ -- 



J.2 QQ2 J. Qj. 



while 



d {rvr 
dr 



-dO. 



A solution of the inhomogeneous equation ( l30l) is 

Vr{r,0,t) 



AMj+pr) BKi(j-(3r) 
Vq H H + 



QqUT 



Qour 



iC 



cosOe 



—iu)t 



(30) 



(31) 



(32) 



with Bessel function Ji of first order and second kind, modified Bessel function Ki 
of first order and second kind, and 



1 



r 



x/2 



(1±^) , /3 



(33) 



The constants A, B and C are still to be determined. The general solution of (|30l) is 
the sum of the general solution of the homogeneous equation plus a special solution of 
the inhomogeneous equation. If, however, the boundary conditions can be satisfied by 
the inhomogeneous solution as given above we need not care about the homogeneous 
solution any more. 

The Bessel function of the second kind Ji{j^ (3r) with complex argument 



j^l^r = ^ (1 + i) J—r 
V2 V 



(34) 



6 



diverges to infinity as r ^ oo. Since the velocity has to be finite at infinity we must 
set A = 0. At the cyhnder surface we have the no-shp boundary condition 



Vr{a, 9, t) ■■ 
and we therefore obtain 



BKAj-f5a) iC 



Qoua 



COS! 



-iujt 







B 



QqO^UJVq + iC 



aKi{j-(3a) 

Substituting the constant B into fl32l) we get 



Vr{r,9,t) =\vo + 



iC 



-[vo + 



iC \ K.irpr) 



Qqujo? J Ki{j (3a) 



cos( 



-ILOt 



while fl3ip and a little algebra provide us with 



veir,e,t) 



iC 



+ 



Vo + 



iC \ K.ij-pr) 



Qo^o? ) Ki{j-(3a) 



(35) 



(36) 



(37) 



(38) 



+a/3j~ 



Vo 



Qo^Q^ ) Ki{j (3a) 



smi 



-iu)t 



There is one last constant C that still has to be determined. We solve the no-slip 
boundary condition vg{a,9,t) = at the cylinder surface for C and find 



C = iQoa^ujvo ( 1 + 
As an abbreviation we introduce 

/(r) : = 

and with 

/? 

we finally obtain the full solution for r > a, 



2j+ K,{r/3a) 
Pa MrPa) 

2j+K,{j-(3r) 
MrPa) 



Vrir,9,t) = t;o cos(^)e-*"* 
ve{r,9,t) = vosm{9)e-''^' 
p{r,9,t) = ivogou cos{9)e 



Koij-Pr) f{r) afja) 



-iujt 



"r2+Vo(j-/3a) ' Pr 

'r + — + ''^^''^ 
r 



Pr 



(39) 
(40) 
(41) 

(42) 
(43) 
(44) 



Figure [2] shows the pressure distribution around the cylinder. Figures [3l HJ and [5] 
contain plots of the fiuid velocity. The fiuid velocity always reaches its undisturbed 
value Vq for r — oo, as required by the boundary conditions at infinity. For high fre- 
quencies, the fiuid velocity reaches its undisturbed value for much shorter distances 
than for lower frequencies. 
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Figure 2: Pressure distribution ^[p{r, 6,t = 0)] around a cylinder (red) with radius 
a = 1 fim for frequency / = UJ/27T = 10 Hz in arbitrary units. The pressure decreases 
steeply within a few cylinder diameters. We assume dry air in all figures. 
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r (mm) 



Figure 3: Amplitude of radial relative velocity \vr{r,0 = 0)|/fo of (H2|) in depen- 
dence upon the distance r to a cylinder with radius a for frequency / = 00/211 = 
10 Hz. The relative velocity approaches 1 for r — cx), as required by the boundary 
conditions. The thicker the cylinder, the farther the influence of the cylinder on the 
velocity field. 
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Figure 4: Amplitude of radial relative velocity \vr{r,0 = 0)|/fo of P2|) in depen- 
dence upon the distance r to a cylinder with radius a = 1 fim for different frequencies. 
For high frequencies, the fluid velocity reaches its undisturbed value at much shorter 
distances than for lower frequencies. 
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Figure 5: Amplitude of relative velocity \vg{r,9 = 7r/2)\/vo of (H3!) in dependence 
upon the distance r to a cylinder with radius a = 1 fim for different frequencies. As 
the fluid must flow around the cylinder, relative velocities larger than 1 occur. For 
high frequencies, the fluid velocity reaches its undisturbed value at much shorter 
distances than for lower frequencies, as in Fig. HI 



11 



3 Force on a Cylinder 



The force / acting on a surface element is given by stress tensor 11 operating on the 
normal vector n of the surface, 



f = Un. 



In polar coordinates the stress tensor is given by (ILambl . 119321 . Art. 328a) 



n 



-p + 2iiodrVr /io ildeVr + drVe - \ve) 



(45) 



(46) 



The normal vector on the cylinder surface reduced to a circle is 



n 



cos 6 
sin 6 



(47) 



Because of ( H5l) and ( l46i) . the force fx per area in x-direction is given by 

fx = {—P + 2flQdrVr) cos 6 — /Iq ( -dgVr + drVg Vg] Sm6 . 



(4J 



The total force Fy in y-direction is zero because of symmetry. The total force in 
x-direction per cylinder length is 



2lT 



-p + 2/io(9rf r) COS — { -dgVr + OrVg Vg ] siu 9 

' ly ry ' 



a dO . 



Substituting of the solution fl42] - |44j) and integration gives 



'i27rgoVoUJa 



a + 



(49) 



Figure [6] shows the force acting on a cylinder in dependence upon the fluid fre- 
quency. For low enough frequencies and thin cylinders forces that are independent 
of viscosity, so-called "buoyancy" forces, 



-i2TTgoVoU}a e 



2 „—iuit 



(50) 



can be neglected as compared to viscous forces. When calculating the force on an os- 
cillating cylinder with angular f requency u additional buoyancy forces proportional 
to qoVquo^ occur (jPantonl . |2005| ). If these are negligible, however, we need only con- 
sider the relative velocity Vq between cylinder and fluid even in case of an oscillating 
cylinder so as to find 



Fx^-t27rgoVou;^e-^-' 



(51) 
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Figure 6: The amplitude of tlie force |-Fx|/^^o per cylinder length relative to fluid 
velocity vq acting upon a cylinder that stands still; cf. (H9l) . The force is plotted for 
different cylinder radii a in dependence upon the frequency /. For a thin cylinder 
(a = 1 /im) the force is more or less independent of the frequency in the range 100- 
1000 Hz. For thicker cylinders (a = 10 /xm) the force increases with frequency. The 



"buoyancy" 
forces. 



force Fp (green graph) of fl50|) is negligible in comparison to viscous 
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4 Discussion 



For a very long cylinder in the low Reynolds number regime time-dependent incom- 
pressible Stokes equations hold. For harmonic flow, these can be solved analytically 
giving the flow and pressure fleld around a cylinder. However, the approximation of 
an inflnitely long cylinder is only fulflUed when the stream velocity converges to its 
undisturbed constant value within a short range, so that the cylinder must be much 
longer than the distance where the velocity has reached about 90 % of the undis- 
turbed value. This precondition is not fulflUed for very low frequencies as shown in 
Figs. H and U 

The velocity fleld and viscous forces can be calculated analytically using the linear 
time-dependent Stokes equations as an approximation. In contrast to Stokes (1851), 
who computed the velocity fleld in the direct neighborhood of a cylinder, and, hence, 
the force for small radius a, we provide the full velocity field for all r > a. 
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